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Abstract. We report the first observations of an astro- 
nomical object using a superconducting tunnel junction 
(STJ) device, a pixel detector with intrinsic energy resolu- 
tion in the optical wavelength range. The Crab pulsar was 
observed using a 6 x 6 array of Tantalum STJs at the 4.2- 
m William Herschcl Telescope on La Palma. Each array 
element provides photon counting capability recording up 
to ~ 10'^ photons pix^^ s^^ with an arrival time accuracy 
of about 5 fis, and providing a wavelength resolution of 
about 100 nm. We derive a photometrically resolved light 
curve which, however, shows no significant colour index 
variations with pulsar phase. 
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1. Introduction 

The p ossibility of intrinsic determination of individual 
photoE energies in the optical range was first reported by 



but is restricted in the present system to about 300- 
700 nm, as a result of the atmosphere 300 nm) and 
the optical elements required for the suppression of in- 
frared photons (~ 700 nm). The detector quantum effi- 
ciency is around 70% across this wavelength range, lim- 
ited by the device/substrate geometry rather than by the 
intrinsic detector response. Count rate limits are about 
10'^ photons s~^, determined by the output stage electron- 
ics, although the device relaxation time is much faster, 
being below ~ 10 /is for the present device. The current 
wavelength resolution, ^ 100 nm at 500 nm, is driven by 
system electronics and residual thermal background (IR) 
radiation, although the intrinsic response of the Ta junc- 
tions is some factor of 5 better than the present perfor- 
mance. 

PSR B0531+21 in the Crab Nebula was first observed 



as an optical pulsar by Cocke et al. (1969) , and provides an 
excellent target for verification of the system's astronomi- 



cal performance. Along with PSR B0833-45 in Vela ( Wal- 



|Perryman et al. (1993), who proposed the application of 



lace et al. 1977|), PSR B0540-69 in the LMC ( [Middleditch 
fc Pennypacker 1985| ) and more recently PSR B0656+14 



STJ te chnology to optical photon counting. Incident pho- 



' tons bi eak Cooper pairs responsible for the superconduct- 



ing state. Since the energy gap between the ground state 
and excited state is only a few meV (rather than ~ 1 eV in 
the case of semiconductors), each individual photon ere 
ates a 



(Shearer et al. 1997) and possibly Geminga (Shearer et al 
1998) it remains one of the few pulsars observed to emit 
pulsed optical radiation. 

While the ~ 33 ms period pulsar has been extensively 



studied at all wavelengths including the optical (e.g. Per- 



gle optical photon counting with energy resolution were 



reported by Peacock et al. (1996) using STJs, and fur 



arge number of free electrons, in proportion to the cival et al. 1993|, Eikenberry et al. 1996, 


Nasuti et al. 1996, 


energy. The first experiments demonstrating sin- Eikenberry & Fazio 1997 




GuU et al. 1998|, Martin et al. 



ther developments have been described by Peacock et al 



(1997). Similar results using superconducting transition- 
edge sensors ( TES) as microcalorim eters have recently 
been reported ( Cabrera et al. 199j), including first ob- 



1998) it continues to be important in providing new in- 
sights into the nature of the pulsar emission mechanism: 
the pulse profile shape, the separation of the primary and 
secondary emission peaks by 0.4 in phase and recent re- 
sults on the energy dependence of the pulse shape over 



(1998) 



the infrared to ultraviolet range (Eikenberry et al. 1996 



servations of the Crab pulsar by Romani et al. 

Technology development within the Astrophysics Divi- 
sion at ESA is ultimately aiming at large arrays of low Tc 
superconductors capable of 10 A intrinsic energy resolu- 
tion at high count rates. A 6 x 6 array of 25 x 25 fim^ Tan- 
talum junc tions has been deve loped as a first astronomical 
prototype ( Rando et al. 199^ ). The wavelength response 
is intrinsically very broad (from < 300 nm to > 1000 nm) 



Eikenberry & Fazio 1997) provide a challenge to theo- 
retical models in which 7-rays created through curvature 
radiation interact with the pulsar magnetosphere to pro- 
duce the X-ray, ultraviolet, optical and infrared pulsations 
through a variety of energy-loss mechanisms. A photon 
counting detector with energy resolution in the optical of- 
fers an important possibility to examine further the energy 
dependence as a function of pulse phase. 
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2. Observations 

Our prototype 6x6 Tantalum STJ array covering an area 
of ~ 4 X 4 arcsec^ was operated at the Nasmyth focus of 
the William Herschel Telescope on La Palma in February 
1999. Photon arrival time information was recorded with 
an accuracy of about ±5 /iS with respect to GPS tim- 
ing signals; while the latter is specified to remain within 
1 lis of UTC, typical standard deviations are much less 
( Kusters 199(]| ). Observations were made on 4-6 Feb, al- 
though modest seeing (> 2 arcsec), especially poor on 
the first two nights, and a significant number of unsta- 
ble junctions, meant that total intensities could not be 
determined reliably. Our present analysis is restricted to 
the signal extracted from just 6 pixels (corresponding to 
an indeterminate but small (~ 0.1) fraction of the overall 
PSF), and a consideration of the resulting pulse profile 
and its energy dependence. 

Data from the STJ are archived in FITS format with 
the photon records defined by their arrival time, x, y pixel 
coordinate, and energy channel in the range 0-255. Chan- 
nels 50-100 cover AA ^ 610 — 310 nm, with A(nm) ^ 
1238.5/(to X A^ch + c), where A'ch is the channel number, 
m ~ 0.04, and c ~ 0.03. Energy calibration was performed 
using an internal calibration source before and after the 
target observations, and verified using narrow-band filter 
observations of a standard star. 

Photon arrival times were translated to the solar sys- 
tem barycentre using the JPL DE200 ephemeris, taking 
into account gravitational propagation delay. Our refer- 
ence timing ephemeris for the Crab pulsar used the 15 Feb 
1999 (MJD = 51224) values oi u = 29.856 514436 4 Hz 
and V = -374 886.90 x 10~^^ s"^ taken from the radio 



ephemeris of Lyne et al. (1999). Consistent periods were 



obtained, with a precision of typically 5 x 10 s, from a 
period search of the timing data. 

3. Results and Discussion 

Fig. shows the light curve from the 50 min of data ob- 
tained on 6 Feb, acquired with a time resolution of 5 /iS, 
and folded into 128 phase bins (~ 250 /iS per bin), with 
an arbitrary origin of zero phase, and without background 
subtraction (impractical due to the combination of poor 
seeing, small field, and fraction of unstable junctions; the 
overall system response is undetermined for similar rea- 
sons, combined with uncalibrated losses at the derotator 
entrance aperture). Examination of the light curve, includ- 
ing the peaks, at finer temporal resolution down to about 
30 /iS per phase bin, reveals no significant sub-structure 
persisting over the observation interval. 

Fig. ^ shows the data divided into two separated en- 
ergy channels, corresponding to Ei = 310-410 nm and 
E2 = 500-610 nm (energy channels 77-99 and 50-61 re- 
spectively). The profiles are normalised to the same rel- 
ative intensities in the peaks, and are displaced verti- 
cally. For this choice of energy bins, the ratio of photons 
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Phase 

Fig. 1. Pulse profile for the Crab pulsar from the 6 Feb 
1999 data over the range 310-610 mn (128 phase bins). 
Photon statistical error bars are included. During the 
10 min interval of best seeing, the total counts summed 
over the 6 active pixels corresponds to about 1000 photons 
s^^, roughly half of which originate from the unpulsed 
component (including sky) and half from the pulsed com- 
ponent. 



is roughly E2:Ei = 4:1. The resulting colour ratio, con- 
structed as E2-E1/E2-I-E1 and folded at the pulsar period, 
is shown in the lower part of the figure. 

Wavelength -dependent variations acr oss the peaks has 

based on obser- 



been noted by Eikcnbcrry ct al. (1996) 



vations spanning the UV to the infrared K band, and by 
^andoval et al. 1998 . Not suprisingly given the form of 
these reported variations, no significant variations with 
pulsar phase are evident in our data. From the nebular 



emission-line spectrum reported by Davidson (1979) the 
[0 11] doublet (A 3726, 3729) lies within our extracted blue 
passband, while the red passband selected comprises none 
of the strong emission lines in the red part of the spectrum 
([Nil] 6584, Ha 6563, [Nil] 6584, [Sin] 6717, 6731), and 
only marginally collects photons from the [O ill] doublet 
(4959, 5007). Our a posteriori choice of energy channels 
for this ratio provides sufficient wavelength separation to 
avoid 'contamination' of each bin given the low energy res- 
olution of the device. Further discrimination of the nebu- 
lar contributions is limited by our presently modest energy 
resolution, which is also insufficient to confirm the reality 
of the b road absorption feat ure near 5920 A, so far noted 
only by [Nasuti et al. (1996) . 
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Fig. 2. Normalised 32-phase bin pulse profiles for the 
6 Feb data (i.e. energy cuts of Fig. 1), divided into the 
El = 310-410 nm (upper curve) and E2 = 500-610 nm 
(middle curve) energy channels, with Acff ~ 378 and 
553 nm respectively. The upper curve has been shifted 
vertically by 1.0. The resulting colour ratio versus phase, 
constructed as E2-E1/E2+E1 (Ei = 310-410 nm and 
E2 = 500-610 nm), is shown in the lower panel; the or- 
dinate is the integrated photon content of each of the 
summed energy channels (arbitrary scale). 



These observational results can be compared with pre- 
dictions from models in which the observed flux in a given 
phase interval is a combination of emission from effec- 
tively disjointed physical regions, in which the sections 
contributing to the emission in a giv en phase interval de- 
pend on the viewin g geometry (e.g. Cheng et al. 1986a , 



Cheng et al. 1986b ). This mixing of physical regions acts 



to average the total emission, and predicts that observed 
properties such as the energy ratio remain constant, or 
change modes tly but rapidly close to the emission peak 
(iRomani 1996| ). Higher S/N ST J data will be required to 
probe the small, rapid spectral index variations across the 



pulse peaks reported by Sandoval et al. (1998 



As the first application in astronomy of a supercon- 
ducting tunnel junction detector capable of providing in- 
trinsic wavelength resolution in the optical, the results are 
a modest indication of the technology's capabilities for the 
future. Significant improvements in performances, and in 
particular in the wavelength resolution, are expected from 
the use of lower critical temperature superconductors in 
the future. 
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